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Abstract 

We analyse the recent precision measurements of the lepton-hadron deep inelastic scat- 
tering at CERN and SLAC to extract model independent constraints among the nucleon 
matrix elements of the twist-4 operators. We also study a parameterization of these matrix 
elements and point out the possibility that the matrix elements of the quark-gluon mixed 
operator has a negative value of the order of —(400 ± 100 MeV) 2 at 5 GeV 2 renormalization 
scale. 



1 Introduction 



Measurements of the lepton-hadron deep inelastic scattering(DIS) remain to be the corner- 
stone of various QCD tests, ranging from precise determination of Aqcd to the knowledge 
of the structure functions necessary to calculate cross sections for hard scattering processes. 
(For the recent review, see Jl|.) 

Recent precision DIS data at CERN || ^ |J and at SLAC || also provide us with a 
fruitful byproduct, i.e. the estimate of higher twist effects in the spin averaged structure 
functions (F2 and Fl). The twist-4 part of these structure functions is defined through 

F 2 ,l(x,Q 2 ) =F 2 ^ 2 (x,g 2 ) + ^F 2 ^ 4 (x,Q 2 ), (1) 



where the target mass corrections H] are taken into account in the twist-2 part F, 



T=2 

2,L ■ 



In terms of the operator product expansion (OPE), the four-quark operators (^r Ml \E'\i/r M2 \I/) 
and the quark-gluon mixed operator (^f{D fll ,*F fl2a }r a ^) contribute to F^j 4 ffi, EH, [9| [10|] . 
The twist-2 part F^ 2 is known to give a parton distribution (i.e. the single particle property 
of quarks and gluons in the nucleon), while the matrix elements of the twist-4 operators are 
the measure of the correlation of quarks and gluons in the nucleon. 

Such new information has wide applications in QCD: first of all it gives a detailed knowl- 
edge of the nucleon structure and gives a stringent test of the various models of the nucleon. 
Secondly, these twist-4 matrix elements are useful to analyse the higher twist effects in other 
high energy processes such as the neutrino induced reaction and the Drell-Yan processes 
[|i~T| . Thirdly, the twist-4 matrix elements are essential to study the propagation of hadrons 
in nuclear medium as is shown in the framework of the QCD sum rules [112 . 



At present, an unambiguous determination of the magnitude of the twist-4 matrix ele- 
ments is not available. However, the recent NMC data 0, |3| together with the SLAC [0] 
and BCDMS || data give us a useful constraint among the twist-4 matrix elements. In this 
paper, we will first examine such constraints in a model independent way. Then, we will 
introduce a parameterization to satisfy the constraints and point out that the quark-gluon 
mixed operator at 5 GeV 2 scale has a sizable nucleon matrix element. 
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2 Operator Product Expansion 



The spin-2, twist-4 contribution to the spin-averaged forward amplitude of the electromag- 
netic current can be written as || |[ 



where the polarization tensors are defined as e^ v = g^ u — q^qu/q 2 and d^ v = —p^p u q 2 /(p-q) 2 + 
(PfiQu +Pv%)/p " Q — g^v with Q 2 =—q 2 . (p^ 1 is a 4-momentum of the nucleon with p 2 = Mjf.) 
A 1,2,9 are the spin-averaged nucleon matrix elements of the spin-2, twist-4 operators: 

(0%) = ( PaP p - \M 2 N g aP )A\ (3) 

with 

0% = g 2 (qi a Q 2 t a q)(q lp t a q), 

0% = tg(q{D a ,*F Pll }ri 5 Q 2 q). (4) 

Here, the operators are assumed to be symmetric and traceless with respect to the Lorentz 
indices: O a p —> \{O a p + Op a ) — \g a pO^. Q is the flavor SU(2) charge matrix and t a are 
the generators of the color SU(3) normalized to tr(t a ) 2 =l/2. F a/3 = F^t a , and the dual field 
strength is defined as *F af3 = e^sF 15 with e i23 = 1- Here we have neglected the twist- 
4 operators proportional to the current quark masses. A typical diagram which generates 
O 1 is given in Fig. 1(a), and that for O 9 , O 2 is given in Fig. 1(b). If one writes eq. (0) as 
T = 2M/x 2 Q 2 , twist-4 matrix elements and the twist-4 structure functions are related as 

M 2 , L (Q 2 )= f dxF^\x,Q 2 ). (5) 
Jo 



3 Experimental data 

Structure Function ^(x) 
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The experimental data of F 2 4 (x) have been analyzed by introducing the following un- 
known function C(x) 

Fr 4 = C(x)F 2 LT (x,Q 2 ), (6) 

where F 2 LT (x, Q 2 ) denotes the leading- twist structure function with the target mass correc- 
tion §. 

C{x) has been extracted for the hydrogen and deuterium target in ref. [F3J by using 



the BCDMS data and the SLAC data taken in the kinematic region 0.07 < x < 0.75 and 
0.5 < Q 2 < 260 GeV 2 . We have carried out y 2 fitting of the proton data C p (x) (given in 
Table 2 of @) by 



C(x) = a + a\x + ci2X + a 3 x + a 4 x , (7) 

and we get a = —0.28, a x = 3.45, a 2 = —17.13, a 3 = 31.64, and a 4 = —14.95. 

One can also extract C n (x) by combining hydrogen and deuterium data in [131]. The 
result, however, has large error bars. On the other hand, the NMC group recently published 
better statistics data for C p (x) — C n (x) (but not for C p (x) and C n (x) separately) which is 
a combination of NMC, SLAC and BCDMS data |J. The NMC group analyzed the ratio 
F£/FP in the kinematic range 0.07 < x < 0.75 and 0.8 < Q 2 < 75 GeV 2 . This ratio is 
independent of the spectrometer acceptance and normalization and gives a reliable estimate 
of Cp(x) — C n (x) from the following relation, 

Fj = { W (1 Q* } - (8) 



By combining this data with that of the proton in ref. ]T3|] and fitting the resulting values for 



C n (x) with the same polynomial in eq. (0), we obtain the following values for the coefficients; 
ao = —0.28, ai = 3.12, 0,2 = —11.01, 03 = 16.51, and 0,4 = —2.40. We checked that different 
set of fittings fall well within the estimated errors of the following results. 

From our fit of C p (x) and C n (x), the integrated structure function at Q 2 = 5 GeV 2 (which 
is a typical scale where the twist-4 effect is extracted) reads 

" A , 1 / . 1 5 . o 1 



/ FZ= 4 dx = -(A 1 + -A 2 + — A 9 ) 
Jo 2 2 V 8 16 ; 

0.011 ±0.004 GeV 2 (neutron) 



1 C(x)F LT (x)dx - \ °-°° 5 ± °-°° 4 GeV2 (pr ° t0n) (9) 
L;{x)i< 2 [x)ax - < nni1 , nnn/l n ,r2 — \ • lyj 



The errors come from unavailability of C(x) for x > 0.75 and x < 0.07. Here we have used 
the leading order (LO) structure function of Gliick-Reya-Vogt for F 2 iT . At Q 2 ~ 5 
GeV 2 , the difference between the LO and the higher order (HO) distribution functions are 
not significant after the x-integration. 1 



Longitudinal Structure Function Fl(x) 

The higher twist effect in the longitudinal structure function is obtained by the ratio 
between the longitudinal and transverse cross sections R = a L /a T . This ratio is especially 
sensitive to the higher twist contribution because the lowest twist effect to Fl is of order 
a s . Note that only diagrams such as given in Fig. 1(b) contribute to Fl. In this case, the 



twist-4 analysis using the transverse basis provides us with an intuitive picture | 10|| , in which 
the higher twist effects can be interpreted in terms of the intrinsic transverse momentum 
of partons: F£ =4 (x) =4/ cPhTk^flx, fcf.), where f{x,k^) denotes a structure function for 
quarks with the momentum fraction x and the transverse momentum /cy. 

Motivated by this, the SLAC data |J was analysed in ref. |15j by introducing a typical 
scale for the transverse momentum of the parton 

FI=\x,Q 2 ) = 8k 2 F 2 lt (x,Q 2 ). (10) 

By using the leading and next-to-leading order structure function for F 2 LT , the SLAC data 
can be fitted by 

k 2 = 0.03 ± 0.01 GeV 2 , (11) 



in the range 0.2 < x < 0.6 [IS, 16). An indirect experimental justification of eq. ( [TD[ ) is that 
R is independent of targets @, ||. If the twist-4 contribution to Fl were not proportional 
to F 2 , the twist-4 contribution to R would depend on the targets. Using the above fit, we 
obtain (at Q 2 = 5 GeV 2 )0 

1 Although this phenomenological parton distribution function might contain the effect of the power 
corrections, this portion will be 0(1/Q 4 ) and thus irrelevant in the twist-4 part of F% in eq. (Q). 
2 This will be an important guide for our parameterization in section 3. 
3 Here we have again used the LO structure function of ref. VM. 



1 %tfF LT (x)dx - I °-° 35 ± °-° 12 GeV ' (pr ° t0n) (12) 
2 [X)X ~ \ 0.023 ± 0.008 GeV 2 (neutron) ' {LZ) 

As is clear from this expression, the difference between the proton and the neutron comes 

only from the difference in / F^dx. 

4 Constraints on A 1 ,2 ,9 

The experimental data for F£ =4: (eq. (fP2l)) is 2-7 times larger than those for F 2 T=4 (eq. @). 
Since both A 1 and A 2 are the matrix elements of the four-quark operators, their absolute 
values are expected to be similar in magnitude. This together with eqs. ( 12 ) and (Q) suggests 



that A 9 at Q 2 = 5 GeV 2 takes large and negative value to reproduce F^ and Fl simultane- 
ously. We will come back to this point in section 5. 

From eqs. ([l|) and @, we can derive two constraints among A 1 , A 2 and A 9 : 

! / -0.165 ±0.061 GeV 2 (proton) 

+ \ -0.093 ± 0.041 GeV 2 (neutron) 

2 3 / 0.280 ±0.096 GeV 2 (proton) 

2 \ 0.184 ±0.064 GeV 2 (neutron) ' 1 j 

The A 1 — A 9 and A 2 — A 9 relations with error bars are given as the bands in Fig. 2. The 

figure shows that it is hard to find a solution where A 1,2 ' 9 are all consistent with zero, which 

clearly indicates sizable values of the twist-4 matrix elements. We note that as long as A 1 , 

A 2 and A 9 do not take too different values among one another, the typical magnitude of 

them reads 0.1 GeV 2 ~ (300 MeV) 2 and a negative value for A 9 is favored. (We will discuss 

this in detail in section 5.) 

Fig. 2 gives an useful test of the various nucleon models: Any reliable models of the 

nucleon should be able to predict the matrix elements within the bands in Fig. 2. One should 

also note that twist-4 data of F 3 (x), although it is not available now, will be particularly 

useful to obtain further constraints on A 1 ' 2,9 . 

5 Parameterization of the matrix elements 



Although Fig. 2 provides us with a model-independent constraint among the twist-4 matrix 
elements, it does not give any definite numbers for the matrix elements. In this section, we 
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will further introduce a theoretical assumption to estimate the magnitude of A 1,2,9 . 
The Bag Model 

The MIT bag model provides us with the simplest estimate of the twist-4 matrix elements. 
Jaffe and Soldate calculated A 1 and A 2 and found that -F 2 T=4 m the model has an opposite 
sign from the data (see the footnote 15 of the latter reference in ||). Shuryak and Veinstein 
H also discussed that models without correlation between quarks inside the nucleon cannot 
reproduce the data. Let's first generalize the MIT bag model parameterization to see whether 
one can remedy the problem encountered in ||. 

The nucleon expectation values of any operators in eq. (||) can be obtained from the bag 
wave function as follows: 



where \p) is the bag state made of three confined quarks. By using the explicit form of \p) , 
one obtains §, A 1 = (2/3)/ia - (16/9)/ 2 a and A 2 = 2j\b + (16/9)/ 2 c. Here the factors 
related to the color- spin- charge read a = — 16/9(— 4/3), b = — 4/3(— 8/9) and c = 8/9(4/3) 
for the proton (neutron). is related to the spacial wave function of quarks: A simple 
estimate with the bag radius 1 fm gives j\ = 0.0266 x a s and / 2 = 0.0042 x a s , which leads 
to A 1 = -0.018(-0.014) x a s GeV 2 and A 2 = -0.064(-0.037) x a s GeV 2 for the proton 
(neutron). a s is the strong coupling constant and we adopt a s ~ 0.5.0 

The mixed condensate can also be obtained from eq. ( |i"4"|) by assuming abelian electric 
and magnetic fields. The electric field vanishes locally within the bag, while the magnetic 
field together with the quark wave function in the bag gives 




(14) 



f 0.075 x a s GeV 2 (proton) 
\ 0.113 x a s GeV 2 (neutron). 



(15) 



Here Oq has a dominant and positive contribution to A 9 . 
Adding all the contributions we finally obtain 




f -0.027 x a s GeV 2 (proton) 
~ \ -0.015 x a s GeV 2 (neutron), 



4 Here it is not obvious whether one should use a s at Q 2 = 5 GeV 2 or something else. In this paper, we 
follow the argument in [p| to estimate an "effective" value a s ~ 0.5. 
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and 

f 1 P r=4 , _ / -0.022 X a s GeV 2 (proton) 

Jo \ -0.026 x a s GeV 2 (neutron). l ' J 

Comparing these with eqs. (^) and (|12|), one finds that the bag model gives incorrect signs 
although the absolute values are the right order of magnitude. The circle (proton) and the 
cross (neutron) in Fig. 2 denote the prediction of the bag model, which shows that the model 
is inconsistent with the current data. 

One may get opposite signs for A 1 and A 2 by making f 2 comparable to f\. However, for 
any reasonable form of the wave function, / 2 is much smaller than fi and in fact the bag 
model gives the most generous estimate. Diquark models give positive signs for the moments 



\T7\, but they do not fit the x dependence of the structure function 



A parameterization based on flavor structure 

Instead of introducing more sophisticated models of the nucleon, we now discuss a differ- 
ent kind of parameterization motivated by eq. (JTQ|) . Let us first rewrite the matrix elements 
of the operators in eq. (§) by using the charge operator Q = diag.(Q u , Qd), 

^p(n) = Q 2 u Ku(d) + Qd^diu) — (Qu — Qd) 2 K ud /2 , 



p(n) — Va^ufd) T ^d ly d{u) 
2 p{n) = Qu^l(d) + Qd^d(u) 
9 p(n) = QuKu(d) + Qd^d{u) 



where if's are the matrix elements defined by 



2 
M 2 
2*9,- 



K ld = -^- 2 {^uY\u){dY\d)) p . 



p- 



(19) 



Here, = lal5 t a , T 2 = la t a , T + = -±(T + T 3 ) and = uY^u + dTJ is a flavor- singlet 
operator. The neutron matrix elements are obtained from those of the proton by the isospin 
symmetry and we have neglected the strangeness contribution to simplify the analysis. 



Noting that the flavor structure of K d 2 ' 9 and that of K^' 2,9 are governed by the <i-quark 
and the w-quark respectively, we will introduce an ansatz in which the ratio K d ' 2 ' 9 / K^' 2 ' 9 is 
equal to the momentum fraction of the d and u quarks in the nucleon: 



IK 



1,2,3 



x(d(x) + d{x))dx/ I x(u(x) + u{x))dx = (3. 



(20) 



Here u(x), d(x), ■ ■ ■ are the usual twist-2 parton distribution functions. (3 takes a value 0.476 
at Q 2 = 5 GeV 2 . The analogous condition for K d ' g /K 2,9 in eq. ( ffif) is a sufficient condition 
to satisfy eq. (|T0[) , which can be checked by substituting eq. flT8|) into eq. ( ]T2|) and equating 
the charge operators in both sides. Thus essentially it does not bring any new constraints. 
On the other hand, the condition for K\jK\ is purely an ansatz: Although it is plausible 
from the point of view of the flavor-structure of the operator, it needs to be checked by a 
non-perturbative method in QCD. 

With eq. (p0|), we can reduce the number of matrix elements from 6 (A 1,2,9 for the proton 
and the neutron) to 4 (K^' 2,9 , K^ d ) ■ Although we have 4 experimental inputs, we cannot 
determine all of them uniquely since the ratio of the proton and neutron data for F[ =A is 
automatically satisfied in our parameterization. Therefore, we will vary K\ d and solve others 
as functions of K\ d . We will also limit the variation of \K\ d \ in the range between \K\\ and 
\K]\. (In fact, the difference between (K^, K d ) and K\ d is only the flavor structure and QCD 
is flavor-blind, therefore these matrix elements should take the similar values in magnitude.) 
The resulting values in GeV 2 unit are given in Table 1. 



K 1 

^ud 


K 


K 2 

u 


K 9 

u 


K 1 

^ud 


K 


K 2 

u 


K 9 

u 




-0.173 


0.203 


-0.238 


-K 


0.083 


-0.181 


-0.494 


(K d + KD/2 


-0.112 


0.110 


-0.300 


-{K d + KD/2 


0.112 


-0.225 


-0.523 


K 


-0.083 


0.066 


-0.329 




0.173 


-0.318 


-0.585 



Table 1 



Table 1 gives the following constraints on the possible range of A 9 at 5 GeV scale: 



(540 MeV) 2 < A 9 < -(340 MeV) 2 (proton) 
-(440 MeV) 2 < A 9 < -(280 MeV) 2 (neutron), 



(21) 



which favor the region inside the parallelograms in Fig. 2. (Note that the results in the 



present parameterization are always confined inside the bands in Fig. 2 contrary to those of 
the bag model.) 

Our analysis here suggests that: 



1. As we have discussed in section 4, the matrix element of the quark-gluon mixed oper- 
ator A 9 is relatively large compared to the four quark operators at 5 GeV 2 scale. The 
magnitude of the former is about —(300 — 500 MeV) 2 which is consistent with a typical 
hadronic scale. The sign and the magnitude of the matrix elements should be under- 
stood in a microscopic manner (either by lattice QCD or by non-perturbative nucleon 
models). To compare model calculations with our results in a quantitative manner, one 
needs to evolve A 1 ' 2,9 from 5 GeV 2 scale to the typical hadronic scale. This requires 
further knowledge of the anomalous dimensions of the operators in eq. (^). Our result 
here is also relevant to the analysis of the QCD sum rules in the nuclear medium |12|| . 

2. One can show that A 1 and A 2 have opposite signs from Table 1. This causes a relatively 
strong cancellation in F 2 T=4 (a;) providing with a reason for the large difference between 
the data on F 2 (eq. ©) and on F L (eq. (0)). 
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Figure Captions 



Fig. 1: Typical diagrams for the twist-4 contribution to the forward Compton amplitude: 
(a) the four quark contribution and (b) the quark-gluon mixed contribution. 

Fig. 2: The twist-4 matrix elements A 1 and A 2 as a function of A g in the unit of GeV 2 . 
They are evaluated at the renormalization scale fj 2 = 5 GeV 2 . The band in solid line 
(dashed line) is a region allowed by the experimental data for the proton (neutron). The 
circle (cross) is a prediction for the proton (neutron) in the MIT bag model. The region 
inside the parallelogram is allowed in the parameterization based on the flavor structure of 
the twist-4 operators. 
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